Intrinsic optical bistability of thin films of linear molecular aggregates: The 

two-exciton approximation 
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We generalize our recent work on the optical bistability of thin films of molecular aggregates [J. 
Chem. Phys. 127, 164705 (2007)] by accounting for the optical transitions from the one-exciton 
manifold to the two-exciton manifold as well as the exciton-exciton annihilation of the two-exciton 
states via a high-lying molecular vibronic term. We also include the relaxation from the vibronic 
level back to both the one-exciton manifold and the ground state. By selecting the dominant optical 
transitions between the ground state, the one-exciton manifold, and the two-exciton manifold, we 
reduce the problem to four levels, enabling us to describe the nonlinear optical response of the 
film. The one- and two-exciton states are obtained by diagonalizing a Frenkel Hamiltonian with 
an uncorrelated on-site (diagonal) disorder. The optical dynamics is described by means of the 
density matrix equations coupled to the electromagnetic field in the film. We show that the one- 
to-two exciton transitions followed by a fast exciton-exciton annihilation promote the occurrence 
of bistability and reduce the switching intensity. We provide estimates of pertinent parameters for 
actual materials and conclude that the effect can be realized. 

PACS numbers: 42.65.Pc, 71.35.Aa; 78.66.-w 



I. INTRODUCTION 

The phenomenon of optical bistability already has 
more than thirty years of history, going back to the theo- 
retical prediction of McCall 1 in 1974, followed by exper- 
imental demonstration of the effect by Gibbs, McCall, 
and Venkatesan^ in 1976 (see also Refs. [HQ], and [H for 
an overview) . Since then a vast amount of literature has 
been devoted to explore the topic (an extended bibliogra- 
phy can be found in our recent paper, Ref. [6J); controlling 
the flow of light by light itself is of great importance for 
optical technologies, especially on the micro- and nano- 
scale. More recently, new materials such as photonic 
crystals^ surface-plasmon polaritonic crystals j& and ma- 
terials with a negative index of refraction^ have revealed 
bistable behavior. 

In our previous work^ we studied theoretically the 
bistable optical response of a thin film of linear molecular 
J-aggregates. To describe the optical response of a single 
aggregate, we exploited a Frenkel exciton model with an 
uncorrelated on-site energy disorder, taking into account 
only the optically dominant transitions from the ground 
state to the one-exciton manifold, while neglecting the 
one-to-two exciton transitions. Within this picture, an 
aggregate can be viewed as a meso-ensemble of two-level 
localization segments,^ which allows for a description of 
the optical dynamics by means of a 2 x 2-density ma- 
trix. Employing a joint probability distribution of the 
transition energy and the transition dipole moment of 
Frenkel excitons, allowed us to account for the correlated 
fluctuations of these two quantities, obtained from diago- 
nalizing the Frenkel Hamiltonian with disorder. By solv- 
ing the coupled Maxwell-Bloch equations, we calculated 
the phase diagram of possible stationary states of the 
film (stable, bistable) and the input-dependent switch- 



ing time. From the analysis of the spectral distribution 
of the exciton population at the switching point, we re- 
alized that the field inside the film is sufficient to pro- 
duce one-to-two exciton transitions, confirming a similar 
statement raised in Ref. fijl 



W 



1> 






* 21 *31^ ^ 
















1 







FIG. 1: Four-level model of the film's optical response. The 
input field induces transitions between the ground state |0), 
one-exciton manifold |1), and two-exciton manifold |2). The 
population of the latter is transferred with a rate w to a 
vibronic molecular level |3), followed by fast relaxation be- 
tween the vibronic sublevels towards the vibronic ground 
state. Finally, the latter undergoes relaxation to the one- 
exciton and/or ground state with the rates 731 and 730, re- 
spectively. The constants 710 and 721 denote the radiative 
decay of the one- and two-exciton states, respectively. 



The goal of the present paper is to extend the one- 
exciton model* by including two-exciton states and tran- 
sitions between the one- and two-exciton manifolds, re- 
spectively. Furthermore, two excitons spatially located 
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within the same localization domain usually quickly 
annihilate, transferring their energy to an appropri- 
ate resonant monomer vibronic level . 12 i 13 ' 14 i 15 ' 16 i 17 ' 18 i 19 
Hence, the generalized model requires the consideration 
of exciton-exciton annihilation.— We will assume that 
exciton-exciton annihilation prevents the three-exciton 
states from playing a significant role in the response of 
the film. The relevant transitions of the model are de- 
picted in Fig. [T] 

To make the two-exciton model tractable, we will 
select the optically dominant transitions between the 
ground state and the one-exciton manifold (as we did 
in Ref. l2(il ). and also between the one- and two-exciton 
manifolds. Treating the different localization segments 
independently, in combination with the state selection, 
allows one to considerably reduce the set of relevant 
states, namely to four states of a segment: the ground 
state, the optically dominant one- and two-exciton states, 
and a high-lying molecular electronic or vibronic state, 
through which the excitons annihilate. This model has 
been implemented for the first time in Ref. ITU using the 
simplifying assumption that the transition energies and 
transition dipole moments are correlated perfectly. Un- 
like Ref. ITl we will account for the correct joint statis- 
tics of both quantities, similarly to our previous work. 6 
The optical dynamics of a single localization segment is 
described within the framework of a 4 x 4-density ma- 
trix. We derive a steady-state equation for the out- 
put field intensity as a function of the input intensity. 
The field inside the film is calculated taking into ac- 
count the field produced by the aggregate dipoles. We 
find that, counterintuitively, tuning away from the res- 
onance may, depending on the dephasing rate, promote 
bistable behavior. In addition, we show that fast exciton- 
exciton annihilation combined with slow relaxation from 
the high-lying vibronic level enhances the tendency to- 
wards bistablity. The phase diagram of bistability is 
computed and compared with the one-exciton model. We 
address also the realizability of the bistable behavior in 
actual materials. 

The outline of this paper is as follows. In the next 
section we present our model of a single aggregate, con- 
sisting of a Frenkel Hamiltonian with uncorrelated on- 
site energy (diagonal) disorder (Sec. Ill Af . Next, we de- 
scribe the selection of the optically dominant transitions 
in Sec. IIIBI and introduce our model for the exciton- 
exciton annihilation of two-exciton states in Sec. Ill Cl In 
Sec. IIID| we formulate our approach, based on the den- 
sity matrix equations in the 4x4 space of states, as well 
as the Maxwell equation for a thin film of oriented linear 
J-aggregates. Section IIIII deals with the results of our 
numerical analysis of the bistable optical response of the 
film in a multidimensional parameter space. We identify 
conditions that are most favorable for bistable behavior 
of the film. In Sec. IIVI we estimate the driving param- 
eters and the input light flux required for experimental 
realization of bistability for films of pseudoisocyanine J- 
aggregates. Section fVl summarizes the paper. 



II. MODEL 

The geometry of the model system and the assump- 
tions we adopt hereafter are essentially the same as in our 
previous paper £ In short, we aim to study the transmit- 
tivity of an assembly of linear J-aggregates arranged in 
a thin film (with the film thickness L small compared to 
the emission wavelength A' inside the film) and aligned in 
one direction, parallel to the film plane. The aggregates 
in the film are assumed to be decoupled from each other; 
their coupling to the environment is treated through phc- 
nomenological relaxation rates (see Ref. @ for a detailed 
discussion) . 

A. A single aggregate 

We model a single aggregate as a linear chain of N 
three-level monomers. The two lower states are as- 
sumed to form multi-exciton bands, as a result of strong 
dipole-dipole excitation transfer interactions between the 
monomers. To simplify the treatment of the multi- 
exciton states, we restrict ourselves to nearest-neighbor 
interactions. The transition dipole moments between the 
two lower molecular states are considered to align in one 
direction for all monomers. Then the (Frenkel) exciton 
part of the aggregate Hamiltonian reads 

N N-l 

H = J2 tnb ] n b n + bl +1 b n ) , (1) 

n—l n—1 

where 6^ (&„) denotes the creation (annihilation) Pauli 
operator of an excitation at site n. The monomer ex- 
citation energies e„ between the two lower states are 
modeled as uncorrelated Gaussian variables with mean 
eo and standard deviation a. The parameter J repre- 
sents the magnitude of the nearest-neighbor transfer in- 
tegral. We assume that it does not fluctuate. After 
applying the Jordan- Wigner transformation, the multi- 
exciton eigenstates are found as Slater determinants of 
one-exciton states ip vn with different z/ . 21 i 22 ' 23 The multi- 
exciton eigenenergies are given by Ylv=i n v £ vi with e v 
being the one-exciton eigenenenrgies and n v = 0, 1 de- 
pending on whether the vX\\ state is occupied or not. 
Particularly, we will be interested in the one- and two- 
exciton states: 

N 

\v) = <p un \n) , ( 2a ) 

N 

\v\vi)= ^2 (VvwiVvina ~ Vvi^P^m) |nm 2 ) , (2b) 

n±>n 2 

where |n) = frJjO) and \nin 2 ) — frj^bJjjO), and |0) is 
the ground state of the aggregate (with all monomers in 
the ground state) . We also will need the transition dipole 
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moments from the ground state |0) to a one-exciton state 
\v) and from a one-exciton state \ v) to a two-exciton state 
| v\ U2}. In units of the single- molecule transition dipole 
moment, they obtain the dimcnsionlcss form 



N 

£ 

71 = 1 



(3a) 



n 2 >n 1 



(3b) 



where it was assumed that the aggregate is small com- 
pared to an optical wavelength. 
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B. Selecting the dominant exciton transitions 

At low temperatures, exciton states reduce their ex- 
tension from the physical size of the aggregate to much 
smaller segments as a result of the disorder-induced An- 
derson localizatio n 24 ' 25 . We will denote the typical size 
of these segments as N* , often referred to as the number 
of coherently bound molecules or localization length in 
terms of the localization theory. 

For J-aggregates, the optically dominant localized 
states reside in the neighborhood of the bottom of the ex- 
citon band. Some of them resemble s-like atomic states: 
they consist of mainly one peak with no node within the 
localization segment (see Fig. |2Ja) . We will denote the 
subset of such states as S. To find all the s-like states 
from the complete set of wave functions ip un , we used 
the rule proposed in Ref. [2ol | ip vn \<Pvn\ | > Co with 
Co = 0.75. The inequality selects those states that con- 
tain approximately 75% of the density in the main peak. 
We found numerically that for a wide range of the dis- 
order strength a (0.05 J < a < J), the thus selected 
states accumulate on average 73% of the total oscilla- 
tor strength (equal to N). Recall that for a disorder- 
free aggregate, the optically dominant (lowest) exciton 
state contains 81% of the total oscillator strength of the 
one-exciton transitions (see, e.g., Refs. I27I andl28h. Fur- 
thermore, we have shown that the s-like states, used as 
a basis to calculate the linear absorption spectrum, well 
reproduce its peak position and the shape of its red part, 
failing slightly in describing the blue wing, where higher- 
energy exciton states contribute as welli^S From this, we 
conclude that our procedure to select the optically dom- 
inant (s-like) one-exciton states works well. 

Similar to the s-like states, one may also distinguish 
states that resemble atomic p states. They have a well 
defined node within localization segments and occur in 
pairs with s-like states. Each pair forms an sp doublet 
localized on the same chain segment. The levels within 
a doublet undergo quantum level repulsion, with their 
spacing nicely following the one that exists between k = 



FIG. 2: (a) The lowest 12 one-exciton states of a chain of 
length N = 500 for a particular disorder realization at the 
disorder strength a — 0.1J. (b) A subset of s states (black) 
and p 3 states (gray) that mostly contribute to the ground 
state to one-exciton and to the one-to-two exciton transitions. 
The average single molecule transition energy eo was chosen 
as origin of the energy scale. 



1 and k = 2 exciton states in a homogeneous chain of 
size N*2§- From the theory of multi-exciton transitions 
in homogeneous aggregates^ we know that the Slater 
determinant of the k = 1 and k = 2 states forms the 
two-exciton state that predominantly contributes to the 
two-exciton optical response. This gives us a solid ground 
to believe that the s-like one-exciton states and the two- 
exciton states composed of (sp) doublets dominate the 
one-to-two exciton transitions in disordered aggregates 
(see below). 

Usually, well defined (sp) doublets occur below the 
bare exciton band edge at the energy —2 J. These dou- 
blets are responsible for a hidden level structure of the 
Lifshits taili^S For the s-like states located close to or 
above the bare band edge, it is already impossible to 
assign a p-like partner localized on the same segment: 
higher-energy states have more than one node and spread 
over segments of size larger than N* (see Fig. [2). To ob- 
tain all the states that give a major contribution to the 
one-to-two exciton transitions, the following procedure 
has been used. First, we selected all the s-like states, as 
described above. After that, we considered all the two- 
exciton states |sz/) given by Eq. (|2b[) and calculated the 
corresponding transition dipole moments /Lt s „. s . From the 
whole set of /x sl /, s , we took the largest ones denoted by 
Msp s ,s , were the substrict s in p s indicates its relation with 
the state |s). This procedure catches all true sp s doublets 
and assigns a partner to solitary s-like states, which do 
not necessarily look like real p states. In Fig. \2$, we de- 
picted the final set of the doublets selected from the states 
in Fig. ^Bp, according to the above procedure, which con- 
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tribute mostly to the one- and two-exciton transitions. 

The average ratio of the oscillator strength of the thus 
selected transitions \s) — > \sp s ) and |0) — > \s) turned out 
to be approximately 1.4. For a homogeneous chain, this 
ratio equals 1.57 (then \s) = \k = 1) and \sp s ) = \ki — 
I, fa = 2)). The similarity of these numbers gives sup- 
port to our selection procedure. Even stronger support is 
obtained from computing the pump-probe spectrum, us- 
ing our state selection, and comparing the result to that 
of the exact calculations^ The comparison revealed that 
the model spectrum only deviates from the exact one in 
the blue wing of the induced absorption peak, similarly 
to the linear absorption spectra. 

C. Exciton-exciton annihilation 

As was already mentioned in the Introduction, two 
excitons created within the same localization segment 
efficiently annihilate (the intra-segment annihilation in 
terms of Refs. HH and l32h . Thus, the authors of 
Ref. [li] studied experimentally the exciton dynamics in 
J-aggregates of pseudoisocyanine bromide (PIC-Br) at 
low temperature and found a 200 fs component in the 
two-exciton state decay. They attributed this to the an- 
nihilation of two-excitons located within the same chain 
segment of typical size of N* = 20. We adopt this mech- 
anism for \sp s ) states described in the preceding section. 
Note that 200 fs is much shorter than all other popu- 
lation decay times. Other processes, such as radiative 
decay, occur at times of tens-to-hundreds of picoseconds. 

Two excitons located on different localization segments 
can also annihilate (the inter-segment annihilation in 
terms of Refs. [12). This process, however, is much slower 



as compared to the intra-segment channel^ we neglect 
it. The thermally activated diffusion of excitons accel- 
erates the annihilation of excitons created far away from 
each other. We consider this diffusion-limited exciton an- 
nihilation as irrelevant to our problem, because for bista- 
bility to occur we need the majority of s-like states to be 
saturated (also see Sec. IIV|) . 

It is usually assumed that the annihilation occurs via 
transferring the two-exciton energy to a resonant molec- 
ular vibronic level (see, e.g., Ref. 1 121) . which undergoes a 
fast vibration-assisted relaxation to the ground vibronic 
state. The population collected in this state relaxes fur- 
ther to the one-exciton state 1 1) of the segment or to the 
ground state |0) of the aggregate (cf. Fig []}. In this 
way, one or two excitations, respectively, are taken from 
the system. In summary, a four-level model, including 
the ground, one- and two-exciton states, and a molecu- 
lar vibronic level through which the excitons annihilate, 
should be employed to describe the optical response of 
the film in the two-exciton approximation. 

D. Truncated density-matrix-field equations 



Within the four-level model introduced in the preced- 
ing sections, we describe the optical dynamics of a seg- 
ment in terms of a 4 x 4 density matrix p a /3, where the 
indexes a and (3 run from to 3, where |1) = \s) and 
2) = \sp s ). We neglect the off-diagonal matrix elements 
P30jP31; and /O32, assuming a fast vibronic relaxation 
within the molecular level 3. Within the rotating wave 
approximation, the set of equations for the populations 
p aa and for the amplitudes of the relevant off-diagonal 
density matrix elements R a p, (a ^ (3) readsii 



Poo = -^Pw [QR* W + O*i?io] + 710P11 + 730P33 , (4a) 

Pll = -7lOPll + 721P22 + 731P33 + ^M21 (^#21 ~ ^*^2l) ~ ^PlO (ftR*i ~ ^* R w) , (4b) 

P22 = - (721 + w) P22 ~ jM21 (^#21 - M*R2l) , (4c) 

P33 = -73P33 + Wp 2 2 , (4d) 

-Rio = - («Aio + Tio) Rio - Mio^ {poo - pii) + 7^21^*^20 , (4e) 

^21 = - ^A 2i + T 2 i + ^W^j i?21 - P2lQ {Pll ~ P22) ~ 7^MlO^*-R20 , (4f) 

-R20 = - ^Aio + iA 2 i + r 20 + ^w^j R 20 + ^21^-Rio - ^ipi riR2i ■ (4g) 

I 



Here, 7 10 = 7o|mio| 2 ancl 721 = 7ol/ i 2i| 2 are the ra- the two-exciton state |2), respectively, with 70 denoting 
diative relaxation rates of the one-exciton state |1) and 
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the monomer radiative rate, and fiio and /Z21 being the 
corresponding dimcnsionless transition dipole moments. 
Furthermore, w is the annihilation constant of the two- 
exciton state |2) and 73 = 730 + 731 is the population 
relaxation rate of the vibronic state |3). The constants 
Tin = 710/2 + T and T 2 i = T 20 = 721/2 + T stand for the 
dephasing rates of the corresponding transitions. They 
include a contribution from the population decay as well 
as a pure dephasing part T, which, for the sake of sim- 
plicity, we assume equal for all off-diagonal density ma- 
trix elements and not fluctuating. By A 10 = u>io — u>i 
and A21 = LU21 — &i we denote the detuning between the 
exciton transition frequencies cjio and cj 2 i, and the fre- 
quency Wi of the incoming field. It is worth to notice that 
Eqs. ([!]) automatically conserve the sum of level popula- 
tions: poo + pu + P22 + P33 = 1- 

The quantity fl = d^E/h in Eqs. ((4]) is the amplitude 
E of the field inside the film in frequency units, where do 
is the transition dipole moment of a monomer and H is 
the Planck constant. It obeys the following equation^ 



N s I 

U = Ui + T R-jy \PWRW + P2lR 



21 



(5) 



where fli = doEi/h is the amplitude Ei of the incoming 
field in frequency units, N s is the average number of s- 
like states in an aggregate, and Tr = 2nn do 2 kL/h is 
the superradiant constant, an important parameter of the 
model s 11 ' 31 In this expression, no is the number density 
of monomers in the film, k is the field wave number, and 
L is the film thickness. The angular brackets in Eq. ([5]) 
denote the average over disorder realizations. 

The set of equations forms the basis of our analysis 
of the effects of one-to-two exciton transitions, exciton- 
exciton annihilation from the two-exciton state, and re- 
laxation of the annihilation level back to the onc-exciton 



and ground states on the optical bistable response from 
an ultrathin film of J-aggregates. In the remainder of 
this paper, we will be interested in the dependence of the 
transmitted field intensity \fl\ 2 on the input field inten- 
sity |f^| 2 , following from Eqs. |J1]) and ([5]). 



III. STEADY-STATE ANALYSIS 

A. Bistability equation 

To study the stationary states of the system, we first 
consider the steady-state regime of the film's optical re- 
sponse and set the time derivatives in Eqs. ([U to zero. 
Furthermore, we will mostly focus on the limit of fast 
exciton-exciton annihilation, assuming the annihilation 
constant w to be largest of all relaxation constants and 
also much larger than the magnitude of the field inside 
the film, |f2|. The reason for the latter assumption is 
based on the fact that below the switching threshold, the 
field magnitude \Q\ 2 ~ {jo&*)r where 70 and a* are the 
radiative decay rate of a monomer and the half width at 
half maximum (HWHM) of the linear absorption spec- 
trum, respectively. As 70 <C cr*, the magnitude of the 
field is also much smaller than a*. Above the switching 
threshold, |fi| becomes comparable to a* & The typical 
HWHM of J-aggregates of PIC at low temperatures is on 
the order of a few tens of cm -1 , which in time units cor- 
responds to one picosecond. On the other hand, the time 
scale of exciton-exciton annihilation is 200 femtoseconds 
(see Sec. Ill Cp . This justifies our assumption |0| <C w 
and allows us to neglect i? 20 in steady-state Eqs. (@|, 
because |i? 2 o| ~ |^/(*A 2 i +T 2 o + w/2)\. Within this ap- 
proximation, we are able to derive a closed steady-state 
equation for the Sl-vs-f^ dependence, which reads 



J 



Ns I 2 £10 

N \ ^10 r 2 



r 2 i + w/2 



1 + 7i?T?(^0 r2 ^" A 2 (P00-PU)+M21 > , 

1 10 + ^10 (r 2 i+w/2) + A| x 



On - P22) 



N s I 2 A10 , . 2 Agi + w/2 

[Poo - pu) + P21— — ; ^2 — t^- (Pu - P22) 



7 ^\^r ?0 + A ?0 



(r 2i + w/2) + A^ 



(6) 



r 



The steady-state populations are given byii 



Poo - Pu = 



Pu - P22 



1 + (1 + W703/71073) S- 



21 



1 + 2510 + (1 + W703/71073) ^21 + (3 + w/73) S10S21 

S10 

1 + 25i + (1 + W703/71073) ^21 + (3 + w/73) S 10 S 2 i 



(7a) 
(7b) 
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where 



Si 



2710 



in 



A? 



r? 



S21 — 



w/2 



2(721+^) A 2 21 + (T 21 +w/2Y 



(8a) 



(8b) 



The terms proportional to hi Eq. ([6]) describe the 
effects of the two-exciton state, exciton-exciton annihila- 
tion, and relaxation from the vibronic level back to the 
one-exciton and ground states. Equation ([6]) reduces to 
the one-exciton model considered in our previous paper— 
by setting /j,2i — 0. Similarly to the one exciton model, 
Eq. ([5]) contains a small factor N s /N, absent in the ear- 
lier paper, Ref . EH This smallness, however, is compen- 
sated by the 7V s -scaling of the average in Eq. ©: it is 
proportional to {(j4 a + (4i))/ N s ~ 2N/N S > l. 6 Thus, 
the actual numerical factor in Eq. (J5]) is approximately 2. 
We stress that, unlike previous work,— Eq. ([6]) properly 
accounts for the joint statistics of all transition energies 
and transition dipole moments. 

It is worth to notice that the second term in the first 
square brackets in Eq. ([6]) represents the imaginary part 
of the nonlinear susceptibility, while the one in the second 
square brackets is its real part. Hence, we will will refer to 
these terms as to absorptive and dispersive, respectively, 
following the convention adapted in the standard theory 
of bistability of two-level systems in a cavity^ 

We numerically solved Eq. (J6j> , looking for a range of 
parameters (Tr, a*, T, 731, 730) where the output-input 
dependence becomes S-shaped, the precursor for bistabil- 
ity to occur. In all simulations, we used linear chains of 
N = 500 sites and the radiative constant of a monomer 
70 = 2 x 10~ 5 J (typical for monomers of polymethinc 
dyes). The exciton-exciton annihilation rate was set to 
w = 5OOO70, corresponding to an annihilation time of 200 
fs^ The average single molecule transition energy eo was 
chosen as origin of the energy scale. 10000 localization 
segments were considered in disorder averaging. 

Figure [3] shows the output intensity Tout = \^\ 2 / (lo <J *) 
versus the input intensity I m = | 2 /(7o<7*), varying the 
superradiant constant Tr from small to large values to 
find the threshold for Tr at which bistablity sets in. The 
relaxation constants 730 and 731 from the state |3) were 
taken to be equal to the radiative rate of a monomer, 
70, which is the smallest one in the problem under study. 
The incoming field was tuned to the absorption maxi- 
mum A^Q = eq — uji — 2.02J, which naively speaking is 
expected to give the lowest threshold for bibtability (see 
a discussion of the detuning effects in Sec. IIII Cp . The 
other parameters of the simulations are specified in the 
figure caption. As follows from Fig. [31 for the given set 
of parameters the bistability threshold is T R = la* . 
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FIG. 3: Examples of the output-input characteristics, demon- 
strating the occurrence of S-shaped behavior in the film's op- 
tical response. Simulations were performed for a disorder 
strength a = 0.1J, resulting in an inhomogeneous HWHM 
a* — 0.024J. The incoming field was tuned to the J-band 
maximum, A^' = eo — u>i — 2.02J. The population relax- 
ation rates of the vibronic state |3) were taken equal to the 
monomer decay rate, i.e., 731 = 730 = 70, while the dephasing 
constant V — 5OO70. In the plot, the superradiant constant 
r_R ranges from a* to llcr* in steps of a* (left to right). The 
critical value for bistability to occur is seen to be T C R = 7cr*. 



B. Effects of relaxation from the vibronic level 

From the physical point of view, the most favorable 
conditions for bistability occur in the case of slow re- 
laxation from the vibronic state |3), which is populated 
via a fast energy transfer from the two-exciton state \2) 
(fast exciton-exciton annihilation). Indeed, under these 
conditions, all population can be rapidly transferred to 
the state |3), and, accordingly, the system can be made 
transparent easier as compared to the case of the one- 
exciton model. Clearly, faster relaxation from the state 
3} to the ground state |0) will deteriorate the condition 
for the occurrence of bistability, while slower relaxation 
improves the situation. Figure 3] demonstrates this. 



C. Effects of detuning 



As we mentioned in Sec. IIII K[ a naive viewpoint is that 
tuning of the incoming field to the absorption maximum 
is expected to give the lowest threshold for bistability. In 
this section, we show that in general this expectation is 
incorrect. 

In Fig.[5]we plotted the results of our simulations of the 
film's optical response as a function of the detuning off- 
resonance, A10, obtained for two values of the dephasing 
constant T. The disorder strength was set to a = 0.1J, 
resulting in an inhomogeneous HWHM a* — 0.024J. 
From these data, one can distinguish two regimes. First, 
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Detuning to the red Detuning to the red 




Detuning to the blue Detuning to the blue 




FIG. 4: Examples of the output-input characteristics, demon- 
strating the effect of the relaxation rates 730 and 731 from 
the vibronic state |3) on the occurrence of bistability. The 
set of parameters used in the simulations are: a — 0.1J, 
A{q = £q — u)i — 2.02J (tuning to the J-band maximum), 
r = 5OO70, and T R = 10a*. 



for a relatively large T = 5OO70 = 0.02J ~ a* [panels 
(a) and (b)] the film's response behaves according to the 
naive reasoning: the output-input characteristic looses 
its S-shaped form upon a deviation of the incoming field 
frequency from the absorption maximum. In contrast, as 
is observed in Figs. [5fc) and (d), for V = 2O70 <C a*, 
when the absorption width is dominated by inhomoge- 
neous broadening a*, tuning away from the resonance 
favors bistability. 



We note that similar behavior has been found for as- 
semblies of inhomogeneously broadened two-level emit- 
ters placed in a cavity ) 2 i 33 i 34 where it was suggested 
that this counterintuitive frequency dependence results 
from the interplay of absorptive and dispersive contri- 
butions to the nonlinear susceptibility We believe that 
our model exhibits the same spectral behavior because 
only the ground state to one-exciton transitions lead to 
spectral sensitivity. The one-to-two exciton transitions 
and the relaxation from the molecular vibronic level do 
not: the former because of the fast exciton annihilation, 
which washes out all spectral details, and the latter be- 
cause it occurs from a relaxed state. Thus, all spectral 
features of the two-exciton model of the film's bistability 
are driven by the ground state to one-exciton transitions. 
In other words, the one-exciton (two-level) model consid- 
ered in our previous paper— is relevant for explaining the 
observed spectral behavior. In this case, the bistability 



FIG. 5: Examples of the output-input characteristics, demon- 
strating the combined effect of dephasing, T, and detuning 
off-resonance, A10, on the occurrence of bistability. In the 
simulations, the following set of parameters were used: a dis- 
order strength a = 0.1J (HWHM a* = 0.024 J), the exciton- 
exciton annihilation rate w = 5OOO70, the decay rates of the 
intermediate vibronic level 731 = 730 = 70, and the superra- 
diant constant Tr = 10a*. Panels (a) and (b) represent the 
results obtained for T = 5OO70 ~ o* when changing A 10 from 
the absorption maximum at A^q = en — W\ — 2.02 J to the red 
(a) and to the blue (b) in 20 steps of 0.0025J. The lighter 
curves correspond to a larger A 10. Panels (c) and (d) show 
similar results obtained for F = 2O70 <s£ a*. 



equation ^ is reduced to 



1 a 



1R 



N 



a4 



r 2 + a 2 + |r2| 2 r 10 / 7 o 



N./ 2 



no 



10 



+ A 2 



10 



|o| 2 r 10 / 7o 



101 



(9) 



In our further analysis we show that, indeed, the in- 
terplay of the absorptive and dispersive terms in Eq. ([9]) 
is responsible for the counterintuitive spectral behavior. 
First, let us assume that we are far outside the resonance, 
i.e., |Aio| is large compared to the absorption HWHM, 
whether the homogeneous (T* = (Tin}) or the inhomo- 
geneous one (tr*). Then, the dispersive term drives the 
bistability, because its magnitude decreases as |Aio| -1 
upon increasing A 10, while the absorptive one drops 
faster, proportionally to A^q 2 . The critical superradiant 
constant for the dispersive bistability has been reported 
tober^ = 4[r* + (r* 2 + A 2 ) 1 / 2 ] (see, e.g., Ref.0) which 
is reduced to T C R » 4|Ai | in the limit of |Ai | > T*. On 
the other hand, we found within the one-exciton model 6 
that close to the resonance (|A 10 | <§C er*), where the con- 
tribution of the absorptive term is dominant, T C R scales 
superlinearly with the HWHM, namely as (cr*/r*) a r* 
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with a « 1.7. Similar scaling (T C R — a* 2 /V*) has been 
obtained in Ref. for a collection of inhomogeneously 
broadened two-level systems placed in a cavity. 

The superlinear dependence of T C R for the absorptive 
type of bistability is a key ingredient in understanding 
the counterintuitive Aio behavior of the film's optical 
response. Indeed, let |Aio| ^> o* and Tr = 4|Aio|, i.e., 
we are at the (dispersive) bistability threshold. Now, 
let us go back to the resonance, where bistability is of 
absorptive nature. Choose for the sake of simplicity T R = 
cr* 2 /r* as the critical value. If 4|A 10 | > a* 2 /T*, we are 
still above the (absorptive) bistability threshold, while 
in the opposite case bistable behavior is not possible. 
For a* ~ r*, the line width is almost of homogeneous 
nature, and tuning away from the resonance deteriorates 
the conditions for the occurrence of bistability^ In our 
simulations, this holds for the case of T = 5OO70 = 0.02 J 
and a* = 0.024J [see panels (a) and (b) in Fig. [5]. 

To conclude this section, we note that the detuning 
effect found in our simulations is asymmetric with respect 
to the sign of Aio: the behavior of / out versus I m is 
different for the incoming frequency tuned to the red or 
to the blue from the absorption maximum. We believe 
that this arises from the asymmetry of the absorption 
spectrum. 



D. Phase diagram 

In Fig. [5] we plotted the results of a comparative 
study of phase diagrams of the film's response calculated 
within the one- and two-exciton model under the reso- 



2.02J. Presented is 



nance condition, A^J = eq 
the critical superradiant constant T C R versus the quantity 
W 1/2 = o* + T*, where T* = T + (j w )/2 is the homoge- 
neous width of the one-exciton transition. The last term 
denotes the averaged rate of population relaxation from 
the one-exciton state to the ground state, see Eq. (|4b[) . 
Roughly, W1/2 can be interpreted as the HWHM of the 
absorption spectrum accounting for both inhomogeneous 
and homogeneous broadening (through er* and T* , re- 
spectively). The upper (lower) solid curve in both pan- 
els was obtained for the dephasing constant T = 2O70 
(r = 5OO70) and varying the disorder strength a. For 
a given T, the film is bistable (stable) above (below) 
the corresponding curve. To compare these results with 
those calculated under the assumption that the detuning 
is the only stochastic parameter^ we also plotted the 
vs W1/2 dependence taking all the transition dipole mo- 
ments and relaxation constants equal to their averaged 
values (dotted curves). 

One of the principal conclusions which can be drawn 
from the data in Fig.[6]is that a more efficient dephasing 
helps the occurrence of bistability: all curves calculated 
for r = 2O70 lie above those obtained for T = 5OO70. The 
physics of this behavior is simple: as the threshold for the 
absorptive bistability is T C R = (er/r*) Q r* (see Sec. lITTCl) . 
a smaller T* gives rise to a higher threshold value for 
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FIG. 6: Phase diagram of the bistable optical response of a 
thin film in the (r_R, PFi/2)-space, where W1/2 = <J* + T* , 
with F* = r + (710 )/2, is used as a measure for HWHM of 
the absorption spectrum accounting for contributions of in- 
homogeneous and homogeneous broadening (through a* and 
F*, respectively) to the total width of the J band. The data 
were obtained by solving Eq. (|6]) for the input field tuned to 
the J-band center, A^g' = eo — uii — 2.02J and varying the 
disorder strength a. In both panels, upper and lower curves 
correspond to T = 2O70 and V = 5OO70, respectively. The 
open circles and squares represent the numerical data points, 
whereas the solid lines are a guide to the eye. The solid lines 
themselves represent the Wi/2-dependence of the critical su- 
perradiant constant T^.. Above (below) the curve for a given 
r, the film behaves in a bistable (stable) fashion. For compar- 
ison, we also plotted the phase diagram calculated under the 
assumption that the detuning is the only stochastic parameter 
(dotted curves, cf. Ref. Hlh . 



Tr. Thus, adjusting the dephasing constant T*, we can 
manipulate the film's optical response. This conclusion 
has been drawn already in Ref. [lOj within the simplified 
one-exciton model. 

Another observation is that the magnitude of the crit- 
ical superradiant constant T C R is considerably lower in 
the two-exciton model than in the one-exciton approach. 
This was to be expected from the physical reasoning 
which we presented above: a fast exciton-exciton anni- 
hilation combined with a slow relaxation from the high- 
lying molecular vibronic level favors bistability. With- 
out showing detailed data, we note that also the critical 
switching intensity, i.e., the intensity calculated at the 
bistability threshold, is smaller in the two-exciton model 
compared to the one-exciton model. In both models, it 
also decreases upon increasing the dephasing rate. 

Finally, from comparison between the solid and dot- 
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ted curves in Fig. [6l it appears that, surprisingly, testa- 
bility is favored by the fact that also transition dipolc 
moments and relaxation constants are stochastic vari- 
ables and not only the detuning, as was assumed in the 
simplified model of Ref. 11. At first glance this seems 
counterintuitive. However, inspection of changes in the 
absorption spectrum allows to shed light of on this re- 
sult. We found that upon neglecting the fluctuations, the 
absorption spectrum, first, acquires a shift which intro- 
duces an additional detuning off-resonance. Second, the 
shape of the absorption spectrum gets more asymmetric. 
As the film's response is sensitive to both the detuning 
and asymmetry, the combined effect of these changes pro- 
duces the observed big difference between the two sets of 
calculations. In principle, this discrepancy may be re- 
duced by adjusting the detuning; it is impossible, how- 
ever, to correct for asymmetry. Most importantly, this 
comparison shows that to adequately calculate the film's 
optical response, fluctuations of all variables should be 
taken into account. 



IV. THIN FILM OF PIC: ESTIMATES 

In this section we will analyze low-temperature exper- 
imental data of J-aggregates of pseudo-isocyanine (PIC) 
to shed light on the feasibility of measuring optical bista- 
bility in a thin film of PIC. We will focus, in particular, 
on aggregates of PIC-Br studied experimentally in detail 
in Ref s. IHIH and [bH At low temperatures, the absorp- 
tion spectrum of PIC-Br is dominated by a very narrow 
absorption band (HWHM = 17 cm -1 ) peaked at A = 573 
nm and red shifted relative to the main monomer feature 
(A = 523 nm). For these aggregates, vibration-induced 
intra-band relaxation is strongly suppressed (no visible 
Stokes shift of the fluorescence spectrum with respect 
to the J-band is observed). This favors a long exciton 
lifetime, which is highly desirable from the viewpoint of 
saturation, and thus for optical bistability. The lifetime 
of the exciton states forming the J-band in PIC-Br is 
conventionally assumed to be of radiative nature. For 
temperatures below about 40 K, it has been measured to 
be 70 psi2£ 

Within the one-exciton model studied in our previous 
paper ^ we found that the number density of monomers, 
required for the driving parameter Yr/<j* to exceed the 
bistability threshold, has to obey no > 10 19 cm -3 . Such 
densities can be achieved in thin films prepared by the 
spin-coating methodi 37 ' 38 Within the extended four-level 
model considered in the present paper, the critical ratio 
of Tr/u* may be even lower. Thus, we believe that from 
the viewpoint of monomer density, J-aggregates of PIC 
are promising candidates. 

Another important requirement for candidates, poten- 
tially suitable for bistable devices, is their photostability. 
J-aggregates are known to bleach if they are exposed for 
a long time to powerful irradiation. Therefore, it is use- 
ful to estimate the electromagnetic energy flux through 



the film. For the field slightly below the higher switch- 
ing threshold, the dimensionless intensity inside the film 
obeys J QUt = \tt\ 2 / (j <7*) « 1 (see, e.g., Fig©. Using the 
expression for the monomer spontaneous emission rate 
7o = 327r 3 ^/(3fiA 3 ), we obtain £ 2 ut m 327r 3 ficr*/(3A 3 ). 
The electromagnetic energy flux through the film is de- 
termined by the Poynting vector, whose magnitude is 
given by S on t = aE 2 ut /(47r). Being expressed in the 
number of photons S ou t / {Hujiq) , passing per cm 2 and 
per second through the film, this value corresponds to 
5 x 10 21 photons/ (cm 2 s). As is seen from Fig. [31 above 
the switching threshold the intensity inside the film rises 
by an order of magnitude. Hence, above threshold the 
electromagnetic energy flux reaches a value on the order 
of Sout ~ 5 x 10 22 photons/(cm 2 s). 

Furthermore, the typical time r for the outgoing in- 
tensity / out to reach its stationary value is on the order 
of the population relaxation time, which is 70 ps, ex- 
cept for values of J out slightly above (below) the higher 
(lower) switching threshold, where the relaxation slows 
down. 6 This means that typically, a nanosecond pulse 
is enough to achieve the steady-state regime. Bearing 
in mind the above estimates for S ou t, we obtain the 
corresponding flux for a nanosecond pulse S ou t ~ 10 13 
photons/ (cm 2 ns). On the other hand, for a thin film of 
thickness L = A/(27r) and number density of monomers 
tiq = 10 20 cm~ 3 , the surface density is noA/(27r) « 10 15 
cm~ 2 . Combining these numbers, we conclude that only 
one photon per 20 monomers produces the effect, which 
is well below the bleaching threshold. 38 



V. SUMMARY AND CONCLUDING REMARKS 

We theoretically studied the optical response of an ul- 
trathin film of oriented J-aggregates with the goal to ex- 
amine the effect of two-exciton states and exciton-exciton 
annihilation on the occurrence of bistable behavior. The 
standard Frenkel exciton model was used to describe a 
single aggregate: an open linear chain of monomers cou- 
pled by delocalizing dipole-dipole excitation transfer in- 
teractions, in combination with uncorrelated on-site dis- 
order, which tends to localize the exciton states. 

We considered a single aggregate as a meso-ensemblc 
of exciton localization segments, ascribing to each seg- 
ment a four-level system consisting of the ground state 
(all monomers in the ground state), an s-like one-exciton 
state, a two-exciton state constructed as the antisym- 
metric combination of this s-like state and an associ- 
ated p-like one-exciton state, and a vibronic state of the 
monomer through which the two-exciton states annihi- 
late. To select the s- and p-like states, a new proce- 
dure was employed which correctly accounts for the fluc- 
tuations and correlations of the transition energies and 
transition dipole moments, improving on earlier worksi*i 
The optical dynamics of the localization segment was de- 
scribed within the 4 x 4-density matrix formalism, cou- 
pled to the total electromagnetic field. In the latter, in 
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addition to the incoming field, we accounted for a part 
produced by the aggregate dipolcs. 

We derived a novel steady-state equation for the trans- 
mitted signal and demonstrated that three-valued so- 
lutions to this equation exist in a certain domain of 
the multi-parameter space. Analyzing this equation, we 
found that several conditions promote the occurrence of 
bistable behavior. In particular, a fast exciton-exciton 
annihilation, in combination with a slow relaxation from 
the monomer vibronic state, favors bistablity. In con- 
trast, fast relaxation from the vibronic level to the ground 
state acts against the effect. Additionally, a faster de- 
phasing also works in favor of the occurrence of bistabil- 
ity. 

The interplay of detuning away from the resonance and 
dephasing was found to be counterintuitive. When ho- 
mogeneous broadening of the exciton states (associated 
with the incoherent exciton-phonon scattering) is com- 
parable with the inhomogeneous broadening (resulting 
from the localized nature of the exciton states), the de- 
tuning destroys bistability. Oppositely, at a slower de- 
phasing, the bistability effect is favored by tuning away 
from the resonance. We relate this anomalous behavior 
to an interplay of the absorptive and dispersive parts of 
the nonlinear susceptibility, which jointly contribute to 
the overall effect. 

We found that in general, including the one-to-two- 
exciton transitions promotes bistability. All critical pa- 
rameters, such as the critical superradiant constant, driv- 
ing the bistability, and the critical switching intensity 
are lower than in the one-exciton models In addition, 
bistable behavior is easier to reach if the ratio of the inho- 



mogeneous and homogeneous width is reduced. We also 
found that the stochastic nature of the transition dipole 
moments (the aspect in which our model goes beyond 
Ref. ITlh strongly influences the film's optical response. 

Estimates of parameters of our model for aggregates 
of polymethine dyes at low temperatures indicates that 
a film with a monomer number density on the order of 
10 20 cm -1 and a thickness of \/2ir, achievable with the 
spin coating method^ is sufficient to realize the effect. 
Under these conditions, one photon per 20 monomers 
produces the switching of the film's transmittivity. 

To conclude, we point out that a microcavity filled with 
molecular aggregate a 39 i 40 ' 41 i 42 ' 43 in the strong coupling 
regime of excitons to cavity modes is another promising 
arrangement to realize an all-optical switch. The recent 
observation of optical bistability in planar inorganic mi- 
crocavities^ and the prediction of the effect for hybrid 
organic-inorganic microcavities 4 ^ in the strong coupling 
regime suggest that organic microcavities can exhibit a 
similar behavior. 
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